コード
s1 <- runif(6 * 6) |> matrix(6, 6) |> scale(center = FALSE)
s2 <- runif(4 * 6) |> matrix(4, 6) |> scale(center = FALSE)
# コサイン距離
d <- 1 - proxyC::simil(s1, s2, method = "cosine", use_nan = FALSE)
w1 <- (\() { x <- sqrt(rowSums(s1^2)); x / sum(x) })()
w2 <- (\() { x <- sqrt(rowSums(s2^2)); x / sum(x) })()
transport::wasserstein(w1, w2, p = 1, costm = d, prob = TRUE)[1] 0.1312392